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Abstract. For the model of a linearly driven quantum anharmonic oscillator, the role of 
damping is investigated. We compare the position of the stable points in phase space obtained 
from a classical analysis to the result of a quantum mechanical analysis. The solution of the full 
master equation shows that the stable points behave qualitatively similar to the classical solution 
but with small modifications. Both the quantum effects and additional effects of temperature 
can be described by renormalizing the damping. 



In recent years, driven anharmonic oscillators have generated a large amount of interest 
because of their use as qubit readout devices (i.e., the Josephson bifurcation amplifier) [1, 2j. 
All these devices operate close to or in the quantum regime. There are two possible stable 
states in phase space, which can be distinguished by their respective amplitude and phase. 
Many studies have been performed on the transition between the stable states O HI [5l O [7] and 
further effects like multiphoton resonances [8l [9] and dynamical tunneling [10^ [Ij make driven 
anharmonic oscillators an ideal model to study thermo-dynamics and quantum effects in a non- 
equilibrium system. Most bifurcation amplifiers are operated in the limit of strong damping to 
speed up the classical read-out process of the amplifier. In this paper, we define an equivalent 
to the classical stable state derived from a quantum master equation and discuss its position in 
phase space as a function of damping. 

The driven anharmonic oscillator is described by 

Hs{t) = ^ + - jx' + F{t)x, (1) 

where F{t) = 2Focos(z^t) stands for the external driving force. In the rotating frame and after 
a rotating wave approximation (RWA) the Hamiltonian reads 

H = Aa^a + xa^a{a^a + 'i-) + f{a^ + a). (2) 



Here, we have introduced the detuning between the oscillator's natural frequency and the driving 
frequency A = 57 — i/ , x = —S'jh/{2m'^Q^), f = FQ/{fiVt)yJfi/{2mVt). We then introduce the 




Figure 1. (a) A plot of the quasienergy g in pliase space. Tlie extrema correspond to the two 
stable states, i.e., the high amplitude state (M) and the low amplitude state (m). The saddle 
point corresponds to an unstable state (s). The three kinds of classical orbits are also depicted 
by black closed curves, (b) A cut through the quasienergy potential g ioi P = 0. The quantum 
levels nearest to maximum (M), minimum (m) and saddle point (s) are shown in blue, red 
and green short lines respectively. Shown in long black line is the quantum level on the outer 
torus, which is close to and can be generate with level m. (c) Squared amplitude of the wave 
functions corresponding to quantum levels shown in Fig. 1(b). We use the same colors and labels 
to distinguish them. The two subfigures correspond to different parameters, i.e., wave function 
on the outer torus is degenerate with the low amplitude state (top subfigure) and non-degenerate 
(bottom subfigure). We used A = 0.027 and /3 = 0.0341 for this plot (degeneracy can be created 
by only minor changes of (3). 

position operator Q and momentum operator P in the rotating frame Q = yj\/2{a) + a), P = 
i^\/2{a) — a), where A = —x/^ is the effective Planck constant. Finally, the scaled quasienergy 
Hamiltonian is defined as H = A'^g/x, where 

Here, we have neglected some constant terms. Now all the properties of this system are 
dependent on one single quantity /3 = —Ij'^xl^- The quantum properties are contained in 
the commutation relation [Q, P\ = i\. The transition to the classical regime is equivalent to 
the limit A — )• 0. The dissipative environment is taken into account using the Lindblad master 
equation 

p(r) = p] + r^{{\ + n)V\a\p + nD{a^\p\, (4) 

where the Lindblad operator is defined through T>\A\p = ApA^ — ^{A^ A, p}, n is the Bose- 
Einstein distribution, and r/ is the dimensionless damping strength [G', '12]. 

Let's first discuss the classical dynamics of the driven Duffing oscillator in the weak-damping 
limit. For < /? < 4/27, the Hamilton function exhibits three extrema in phase space, as 
indicated in Fig. 1(a), i.e., a maximum a minimum m and a saddle point s. Without friction, 
there are three kinds of possible periodic orbits. They correspond to different isoenergetic 
sections of the Hamilton function. The orbits with an ear-shape around the maximum M form 
one group. Those circling the minimum m form another group. The third group consists of the 
orbits on the outer torus. One can calculate the explicit expressions of these orbits based on the 
classical equation of motion: dtQ = dpg, dtP = —dqg. In the vicinities of the extrema M and 
m the system is equivalent to an underdamped harmonic oscillator. Therefore, if damping is 
included the orbits nearby will shrink towards the extrema. As a result, they are stable points. 




Figure 2. (a) Positions of the low amplitude state (red solid line), high amplitude state (blue 
solid line) and the unstable saddle point (green dashed line) from the solution of the classical 
equation of motion for finite damping and from the solution of the master equation (circles) 
respectively. The black empty circles represent zero temperature, while the yellow and red solid 
circles correspond to temperatures T = 2.0r2 and T = 3. Oil respectively, (b) The difference 
between the effective damping r/cfj and real quantum simulation damping r/ at zero temperature 
(red solid circles) and T = 0.5f] (black solid squares) respectively, (c) The relationship between 
effective damping and temperature at fixed damping ry = 0.03. Other parameters are A = 0.027 
and /3 = 0.12. 

corresponding to the oscillations with high and low amplitudes respectively. They are separated 
by a phase space barrier associated with the unstable saddle point s. 

In the quantum limit, the classical orbits will be quantized into a series of discrete energy 
levels. They can be calculated by diagonalizing the scaled Hamiltonian ([3]), which results in 
a set of eigenvalues and eigenstates, i.e., g\n) = gn\n). In Fig. 1(b) we show a cut through 
the quasienergy potential g for P = 0. The classically stable states are the extrema of the 
potential. The quantum levels corresponding to the three extrema are shown by colored short 
lines and labelled by M,m and s respectively. The black long line indicates one quantum level 
on the outer torus, which is close to and can be resonant with the level m. Fig. 1(c) shows the 
quantum analogues (i.e., squared amplitude of the wave functions, \^\'^) of the classical orbits 
corresponding to the quantum levels shown in Fig. 1(b). We use the same colors and labels in 
each figures. In principle, the small amplitude state is coupled via tunnelling to states on the 
outer torus. If a state on the outer and inner torus are degenerate, the eigenstates are mixed 
states and this results in a large change of the degenerate wave functions (see the red and black 
curves in Fig. 1(c)). However, one should note that in most of the relevant parameter regime the 
dynamics of the system is dominated by intra- well transitions O [S] . 

As damping increases, the positions of the stable points in phase space will shift. We can 
find the classical solution for the stable points from the classical equation of motion including 
damping: Q = dg/dP — rjQ, P = —dg/dQ — rjP. For finite damping the system is classically 
bistable for ^(i) < /3 < /3^^\ with = 2(1 + V ^ (l _ 3??2)3/2)/27 [6j. In Fig.2(a), we plot 

the possible positions of the low amplitude state (red solid line), high amplitude state (blue solid 
line) and the unstable state (green dashed line) in phase space. 

In order to describe the damping effects in the quantum regime, we turn to the master 
equation (jj]). We can also define a stable state using our master equation. At exactly zero 
temperature, it is possible to get an analytical expression for the stationary distribution in the 
so called P-representation [131 HI] > which is one of several quasi-probability distributions [151 US] • 
With the help of the P-representation, the exact analytical solution for the stationary density 
matrix for n = can be obtained in the photon number representation (see Eq.(15) in Ref.[14j). 
We can diagonalize it and use its eigenstates as our basis. A diagonal representation will not only 
give us a clear physical picture of the dynamics of the system but also simplifies any numerical 



simulation. 

To understand how our solution compares to the solution of the classical equation, we 
use the following procedure: (1) diagonalize the exact solution; (2) identify the eigenstates 
which correspond to the two stable states and calculate the expectation value of position and 
momentum (or amplitude and phase). An example for the results of such a calculation are the 
black empty circles in Fig. 2(a); (3) compare the expectation values of the wave function with 
those expected from the classical equation of motion. Since the states are located on the same 
curve as the classical state we can assign each state a value r/efr that corresponds to the damping 
of the equivalent classical stable point. 

Due to the finite effective Planck constant A used in this work, our results display some 
quantum effects. One is the disappearance of the small-amplitude state already for weaker 
damping than classically predicted. The reason is that near the bifurcation point, the low 
amplitude well is so shallow that the zero point energy exceeds the potential barrier. As a 
result, there is no quantum level in the vicinity of the minimum m. The situation is similar 
for the large-amplitude state if /? is small. Another quantum signature is the renormalization 
of damping, rjcs 7^ rj. We show a comparison of rj, as chosen for the master equation and the 
effective damping rj^s of the equivalent classical stable point in Fig. 2(b). As A becomes smaller 
we find ^ff ~ V %d2i- 

The effect of temperature can be included in an exact numerical solution of the master 
equation. Then we again follow the steps outlined above and find that the eigenstates that 
diagonalize the matrix have average momenta and coordinates that fall on the curve of the 
classically stable states (see yellow and red solid circles in Fig. 2(a)). We find that the 
temperature also introduces a weak renormalization of the damping rj. We plot the difference, 
"HcS — V-i function of rj in Fig.2(b) at a finite temperature (black solid squares in Fig.2(b)). We 
also show the renormalized damping as a function of T for fixed 77 in Fig. 2(c). As one can see, 
there is a small quantum correction for T = 0, and temperature has, as expected, no effect as 
long asT <^VL. As temperature increases we get an additional small shift along the curve of the 
classically stable state. Generally we find that the states that diagonalize the density matrix at 
r = still remain an efficient diagonal basis up to T « O. 
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